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I. INTRODUCTION 



In the standard model of the particle physics (SM), neutrinos are massless particles due 
to the absence of right-handed neutrinos z/r. The simplest way to give masses to three 
neutrinos is to add three Upt similarly to other fermions, which corresponds to six additional 
particles (three and three W^) to the SM. In the Higgs triplet model (HTM) which 
we deal with in this article, a complex SU{2)i triplet scalar with the hypercharge Y = 2 is 
introduced to the SM in order to have neutrino masses. This model can be regarded as one 
of the simplest extension of the SM because the number of new particles is six in this model 
also. 

The triplet Higgs boson field with hypercharge Y = 2 can be parameterized by 

where is the vacuum expectation value (VEV) of the triplet Higgs boson. The constraint 
on the rho parameter, po = 1-00041q;qqq7 at 2a CL (page 137 of 3|]), gives an upper limit 
^^'a/^' ^ 0.01 where v = 246 GeV is the VEV of the doublet Higgs field, which corresponds to 
f ^ < 3 GeV. There is no stringent bound from quark sector on triplet Higgs bosons because 
they do not couple to quarks. The interaction of the Higgs triplet with lepton doublets 
Le = {i^eL, hV = e, p, r) is given by 

Ariplet- Yukawa = L'^ iT2^Lii + H. C. (2) 

The symmetric matrix hu' is coupling strength, Ti{i = 1-3) denote the Pauli matrices, and 
= iyJ^C, iJ^C) with the charge conjugation operator C. 

The coupling hei> has a direct relation to the neutrino mass matrix m«' in the flavor basis 
through as 

= (Umns diag(mi, mse""^!, mae"*'^^) t/nNs) = -^—mu'. (3) 

v2ua ^ v2fA 

The mass eigenvalues rrii are taken to be real positive values. We define Am?^ = mj — m'j 
and refer to the case of Amg^ > (Amg^ < 0) as the normal (inverted) mass ordering. 
Here neutrinos are required to be Majorana particles^, and (fi and (f2 are the Majorana 



^ In general, even if the Higgs triplet exists, neutrinos can be Dirac particles by adding vr also and requiring 
lepton number conservation which results in va = 0. In the HTM we use, neutrino masses are assumed 
to be given solely by va and neutrinos are Majorana particles by definition. 
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BE 



a 



phases [2|, |J] defined in an interval of [0, 27r). Tlie Maki-Nakagawa-Sakata matrix (5|] of the 
neutrino mixing^ is parameterized as 



MNS 



C12C13 



■5l2Cl3 



-S12C23 - Ci2S23Si3e"^ C12C23 - Si2S23Si3e*'^ 



S23C13 



iS 



S12S23 - Ci2C23Si3e' 



-C12S23 - Si2C23Si3e'" C23C13 



(4) 



where s,- 



sin^jj and Cij 



cos 6ij, and 6 is the Dirac phase. The ranges are chosen as 



< Oij < 7r/2 and < 6 < 27r. According to current constraints from neutrino oscillation 
experiments [?, s], o], Q, we use the following values in this article 



Am^^ = 7.6 X 10-^eV^ \Aml^\ 



sin^ 2^12 = 0.87, sin^ 2^23 
sin2 2ei3 < 0.14. 



2.4 X 10"^ eV^ 



(5) 
(6) 
(7) 



The absolute scale of the neutrino mass is constrained by tritium beta decay measurements 



ogical observations as ^ m,, < 0.61 eV (95% CL) 



as < 2.3 eV (95% CL) (llj and by cosmo 
or < 1.3 eV (WMAP only, 95% CL) 

The HTM has seven physical Higgs bosons which are two CP-even neutral bosons 
(lighter) and (heavier), a CP-odd neutral one A^, a pair of singly charged bosons H^, 
and a pair of doubly charged bosons H^^. These Higgs bosons contribute to many lepton 
flavor violating (LEV) processes. Experimental searches for fi eee etc. put upper bounds 
on \h^e\\hee\/'m1[±± etc. (See e.g. [3]), where mH±± is the mass of H^^. The couplings hui 
are, however, not free in the HTM because they relate directly to the neutrino mass matrix 
mw as shown in ([3]). Previous works for dependences of LEV processes on the parameters in 
rriui can be found in . In this article, we consider in detail the correlation of upper 

bounds on \h*jhki\ / m^j^^ from new physics searches and deal with them as lower bounds on 



^ We took the definition of the mixing vi = UaVi according to page 517 of [3| although another definition 
— U^^Vi is used for example, in Q and on page 163 of In latter definition, we need to take 



complex conjugate in the middle equation of ([3]). 
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II. LOWER BOUND ON VAmH±± 



A. Constraint from the Muon Anomalous Magnetic Dipole Moment 



Let us consider first the anomalous magnetic dipole moment (MDM' 



of muon, 



3 



as 



{g — 2)/2. The muon anomalous MDM has been measured very precisely 

a^^P = 11659208.0(6.3) x 10"^°, (8) 

where the number in parentheses shows la uncertainty. On the other hand, the SM predicts 

af^lr] = 11659193.2(5.2) x 10"^°, (9) 
af^le+e-] = 11659177.7(5.1) x 10"^°, (10) 



where the hadronic contributions to ^^^[t] and af^[e~^e ] were calculated ITj 



of hadronic r decay and e^e annihilation to hadrons, respectively (See also 18 



by 



using data 



19 



20 



21 



22|). The deviations of the SM predictions from the experimental result are given by 



Aa^[r] = a;fP - af^lr] = 14.8(8.2) x 10"^°, (11) 
Aa^[e+e-] ^ a°"P - aj^fe+e"] = 30.3(8.1) x 10^^°. (12) 



These values of Aa^[r] and Aa^[e"'"e~] correspond to 1.8cr and 3.7(7 deviations from SM 
predictions, respectively. 

New contributions to at the 1-loop level in the HTM come mainly from H^^ and H^. 
The Yukawa interactions of H^, which are mixtures of doublet and triplet Higgs bosons, 
and of H^^ (= A^^) are written by 

^Sietytkawa = -V2^^==iU^^,h),,l^PJH^ - h,eW'H++ + H.C., (13) 

where Pl = (1 — 7^)/2 and z/j represent mass eigenstates of Major ana neutrinos which satisfy 
conditions Ui = uf. The 1-loop contribution of through the triplet Yukawa interaction^ 



The contribution through m^/v is ignored because it is suppressed by v\lv^ . 
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is calculated as 



m. 



( ^MNS ) ( ^MNS ^) «M 



dt 



K±t' + (1 



K'rr±)t + R 



{rn?) 



MM 



m. 



(14) 
(15) 



and the H^^ contribution is given by 



m. 



V'(/i'^)M^/i^^» / dt 
Jo 



DM 



T-)M 



2e{l-t) 



R 



DM 4^2 



{R' 



H±± 



l)t 



MM 



127r2 v\m]j±± 



(16) 
(17) 



Here we have defined 



m 



2 ' 



U' 



Umns diag(m^, mj, m^) Ul^^g 



(18) 
(19) 



Note that {rn?)ui does not depend on Majorana phases and {m?)u is positive definite. 
Thus a|^™ = a^^ + a^^^ is negative definite though Aa^ is positive. The minus sign of 



the contributions from Higgs triplets has been known 13|, l23|| but it does not seem to be 



taken seriously probably because of confusions about the combination (h^h)^^. Usually the 
combination seems to be written as {h^^Y or (^^)/i/ij for which sign of a^^"^ can be fiipped, 
and then it seems possible to obtain a (finite) constraint on hui. Actually, any value of /i^^/ 
can not fit Aa^[e^e"] (Aa^[r]) at 3.7cr (1.8cr) because of the wrong sign of 

Concerning only on the sign, the 1-loop contribution from can have the right sign 
to explain Aa^ (see 2J] for the case in the type II two Higgs doublet model (2HDM-II)). 
However, the contribution has a suppression with v\lv^ in the HTM because A° does not 



couple with charged leptons at the tree level^. Although in the Barr-Zee type 25| 2-loop 
diagrams the right-sign contribution of can be important in some models like the 2HDM- 



II 



26[ and the minimal supersymmetric standard model (MSSM) [27(| , such a situation does 



The modification of tlic contribution to from /i is also suppressed by v\/v . 
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sin^2e^3 = 0, rn^ = 0, (p^ = 92 = 0. normal, m^++/m^+ = 1 



1000 

5^ 















a with e"*" 
with 
1^ 


3" data 

T data T - 

-» eee — x -> |iee - 
i^ey tie^M'S" 





1 1.5 2 2.5 3 3.5 4 4.5 5 



confidence level [o] 



FIG. 1: Lower bounds on Vi\mjj±± given by constraints on muon anomalous MDM and LFV 
processes as functions of the confidence level. All parameters in the neutrino mass matrix are fixed 
as an example (see the text for the values) for the normal mass ordering. We take mjj±± = mjj± 
for simplicity. 



not happen in the HTM because couphngs of Im(A'')) with quarks and charged leptons 
are also suppressed by va/v. 

The definite sign of is a feature of the simpleness and the predictabihty of the 

HTM. In the MSSM in contrast, the contributions from supersymmetric particles to can 
have the right sign easily by the appropriate choice of the sign of the Higgs mass parameter 



fJ-H 



291]. As the result, the HTM is somehow disfavored by the muon anomalous MDM 



and it results in a strong constraint on the model. This is also the case for other models 
(e.g. the Zee-Babu model 30|) which do not have extra neutral Higgs bosons with sizable 
couplings to charged leptons similarly to the HTM. Of course, the positive Aa^ does not 
seem conclusive yet and it does not mean exclusion of the HTM. The difference between 
Aafj,[e~^e~] and Aa^[r] may indicate existence of new physics in the quark sector which is 
not modified in the HTM. 

Hereafter we take mH±± = rriHi for simplicity. The large splitting of their masses is 
disfavored by the constraint on the p parameter. Once we fix the neutrino mass matrix, 
muon anomalous MDM and LFV processes are interpreted as lower bounds on v/^mfj±±. 
Figure [1] shows the lower bounds with respect to the confidence level in a unit of the standard 
deviation a and they are given by constraints on the muon anomalous MDM with e~^e~ data 
(bold solid red line), the MDM with r data (solid red line), fi — eee (bold dashed green line). 
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/i 67 (dashed green line), r fififi (bold dash-dotted blue line), r Jiee (dash-dotted 
blue line), and the muonium (/i+e") conversion to the anti-muonium (bold dash-dot-dotted 
magenta line). Bounds from r flfifi and r Jiee are important in our analysis among 
six possible r — > Wl" . Formulae of branching ratios of these LFV decays in the HTM and 
their current bounds at 90% CL are 

|2| 



eee 



BR(/i 
BR(/x ^ 67) 
BR(r fififi) 
BR(r /iee) 



< 1.0 X 10 



-12 



3l| 



27a|(m2) 



e/x I 



< 1.2 X 10 



-11 



32| 



MM I 



leCl^-u^m^ii 

I 1 2 1 1 2 

l^^r^l l^eel 

leGl^t'im^ii 



BR(r ^ fiUf^Ur) < 3.2 x 10" 



BR(r ^ /ii^/,z/r) < 2.0 X 10" 



(20) 
(21) 
(22) 
(23) 



where BR(r ^v^Vt) = 17%, a = 1/137 stands for the fine structure constant, and 
Gp = 1-17 X 10^^ GeV~^ denotes the Fermi coupling constant. The effective Lagrangian for 
the muonium conversion is 



C 



MM 



(24) 



The formula of the coupling G^jjj in the HTM and current constraint at 90% CL for that 
are 

2 



\G 



MM I 



Gf 



irrip 



\m 



MM I 



128Gj,vim%^^ 



< (3.0 X lO"'^)^ (34|. 



(25) 



In Fig. [H parameters of the neutrino mass matrix are fixed by ([5]), ([6]), and the following 



values as an example: mi = 0, sin^ 26'i3 = 0, ipi = 



we have {m^)^^i = 1.2 x 10 ^ eV^ 



\m 



fj,e I 



imp 



0. With these values of parameters, 
6.3xlO^^°eV^ 



3.5 X 10^^eV^ \mrf. 



6.4 X 10-^1 eV^ Km2),^|2 
3.6 X 10"^ eV^, and Im^ 



m 



MM I 



5.9 X 



10-9 eV^. Bounds (I20l)-(I23D and ([25]) at 90% CL are translated naively into xa CL bounds 
by multiplying x/1.64 because 90% CL corresponds to 1.6Aa. Below around 1.8cr (3.7cr), the 
muon anomalous MDM Aa^[r] {Aa^[e~^e~]) gives the strongest constraint on the HTM but 
it becomes weaker rapidly than other constraints at higher confidence levels. Hereafter, we 
take Aa^[r] and concentrate ourselves on 2a CL in order to avoid qualitative disagreement 
with Aa^[r] in the HTM. 
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lower bound on ni^++ [eV»GeV] from n -> eee, sin 29^3 = 0, normal 
2 



(a) 



1° 
10, 

10* 
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lower bound on m^++ [eV'GeV] from n - 
2 



eee, sin 29,, = 0, inverted 




m, [eV] 



ms [eV] 



FIG. 2: Contours of the lower bounds on v^mjj±±[eV-GeV] given by /i — > eee. (a) for the normal 
mass ordering, (b) for the inverted mass ordering. 



sin 29^3 = 0, <Pi = 92 = 0- normal, m^^^!m^^ = 1 



sin 29^3 = 0, (pi = 0, 92 = normal, m^^^/m^^ = 1 



(a) 








with X data 








li^ey 









X |iee 




(xe -> |i e - " - " ■ 



(b) 


















2o C.L. with X data 

\i -> eee 

X-> 

X |iee 

|ie ^ |i e - " - " - 
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m, [eV] 



0.25 0.3 



FIG. 3: Lower bounds on v^mH±± for = in the normal mass ordering, (a) for = 0. (b) 
for ip2 = TT. 



B. Constraints in the case of BR(;U — > eee) ^ 

In most of parameter space, the strongest lower bound on VA'mH±± is given by /i — > eee 
as expected naively from the strong constraint on its branching ratio fl20|) . Figures [2](a) and 
(b) show contours of the bounds with ^13 = for the normal and inverted mass ordering, 
respectively. Note that BR(/i eee) does not depend on 6 and for 6*13 = 0. Although the 
bound on VAfnH±± from /i eee is relatively weak for small mi in Fig. [2](a), bounds from 
other LFV processes are weaker than that. It is shown also that ~ makes the bound 
from yU — > eee weak for both of mass orderings. We focus on the case of ipi = in the next 
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|i -> eeevs. ^lee, sin^29^3 = 0, = 0.2eV, normal 

2 
1.5 

0.5 


0.001 0.01 0.1 1 

FIG. 4: In the shaded region, the bound on v/\mu±± from r Jiee is stronger than the one from 
H — > eee. We used mi = 0.2 eV for the normal mass ordering. The shaded region is symmetric 
under a transformation of (<^i, ip2 — t^) ^ {—^i-, —^2 + tt). 

paragraph. In Fig. [2I^a) there is a special point aX (pi = tt and mi = a/ Am2i/ V cos 26^12 — 
4.6 X 10"'^ eV where the bound vanishes because of meg = 0. Such cases of BR(yU — > eee) = 
are discussed in the next subsection. 

In Fig. [3l mi-dependences of bounds from Aa^[r] and LFV processes are presented for 
the normal mass ordering at = where the bound from /i — eee is relatively weak. 
Figures [3](a) and (b) are obtained for (/)2 = and vr, respectively. Other parameters are 
the same values as ones in Fig. [TJ It is seen in Fig. [3](a) that /i eee still gives the most 
stringent bound for ipi = = ^ and mi < 0.3 eV although the bound from the muonium 
conversion gets close to that for large mi. If we accept mi > 0.3 eV, the bound from the 
muonium conversion can be stronger than the bound from /i — > eee. On the other hand, 
Fig. [3](b) shows that the bound from r fiee can be more stringent than the one from 
fi eee for mi > 0.06 eV. This is because a parameter set {9i3,(pi,ip2) = (0, 0, vr) in the 
region of Amfj/m1 ^ 1 gives 

Im^elVeel' ^ ^{Amj^f sin^ 29i2 - 1.6 x 10^^°eV^ (26) 

l^^r^n^ieeP - (27) 

and the large difference between experimental constraints ( l20l) and ( l23l) can be compensated 
for mi > O(0.1)eV. In Fig. H] the shaded region shows values of Majorana phases for 
which the bound from r — fiee becomes more stringent than the one from /z eee at 
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mi = 0.2 eV for the normal mass ordering. The region is symmetric under a transformation 
of {ipi,ip2 — tt) ^ (— v?!, — (/?2 + tt) because of \mu'\ = Although the bound from 

/i — > eee is relatively weak for ipi ^ 0, the bound is still the most stringent one at around 
(/92 = because r Jiee is also suppressed. If we take nonzero ^13 and ignore ^m^^ for 
{^11^2) = (0,vr), eq. is rewritten as 



'^'meel^ ^ 2si3mi (28) 

while eq. (l27j) remains valid. Therefore, the shaded region in Fig. H] at around {ipi,ip2) = 
(0, tt) exists for sin^ 2^13 < 10~^. For the inverted mass ordering, the region where r Jiee 
becomes remarkable is almost same as the one in Figs. [3] and H] because neutrino masses are 
almost degenerated in the reg ion. In such a region, we can also expect a signal of r ^ /2ee 



in future experiments 



35 



with satisfying the current constraint on /i ^ eee. 



C. Constraints in Cases of BR(/x eee) = 



It has been known that the strong constraint from fi — > eee can be evaded in the cases 
of rriefj, = [1^ and = j^. While it is impossible to have = with 6*13 = 0, the 
case of rriee = is possible also for 6*13 = as we mentioned for Fig. [2|^a). Such cancellations 
in the HTM are desired also for experiments 35, 36, 371 to discover some LFV decays 
(yU — > e7 etc.) [6|, [ij] in the future. Figures [5](a)-(d) show results for the case of = in 
the normal mass ordering. Four CP conserving sets of Majorana phases are taken for the 
figures as examples. We use appropriate values of 613 and 6 for rrie^ = 0, which we call as 
"magic values" 6'^'^ and 5™^'^, and explicit formulae of them are shown in Appendix. For 
each cases in Fig. [5l the magic value 5™^'^ is or vr independently of mi. Although {'m^)e^ is 
independent of mi and Majorana phases, the bounds from /i ^ e7 in Fig. [5] are not constant 
with respect to these parameters because O^^'^ depends on them. We see in Figs. [5](a)-(d) 
that the bound from — s> e7 is the strongest one for mi ~ 0, and this is also the case with 
any values of ^pi and if2- For the case of Fig[5t^a), r decays give the most stringent bound 
for 0.01 eV < mi < 0.12eV, and the bound from the muonium conversion becomes the 
strongest one for mi > 0.12eV. In Fig. [5](b), the bound from r — > Jiee is prominent. The 
magic ^13 in Fig. [5]^b) gives sin^ 26'^^'^ ~ 10~^ for mi = 0.2 eV, and then the remarkable 
behavior of the bound from r /2ee can be understood also as a part of the case shown in 
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(Pi = Tt, 92 = 0. normal. m|_|++/m|_|+ = 1 , hgj^=0 



(c) 2o C.L. with T data 

H->eY 

T^HHH 

T -> ^ee 


^^VT^^ 




M/ 




± 





0.002 



0.004 0.006 0.008 

m, [eV] 



0.01 



(p^ = (P2 = 71, normal, mHWmH+ = 1 , hg. 



(d) 2a C.L with T data 

T^tillll 

X -> \iee 

lie ^ |i e 




















'♦.(.* 







0.002 



0.004 0.006 
m, [eV] 



0.008 



0.01 



FIG. 5: Lower bounds on V/\mjj±± for cases of = in the normal mass ordering. The 
value of 6^^^ varies to keep = (See Appendix), (a) for (991,(^2) = (0,0), (J™^'^ = vr. (b) for 
{ipi,^2) = (0,7r),,5-s^ = 0. (c) for (^1,^2) = (vr, 0), = 0. (d) for (^1,^2) = (vr, vr), ^^^^ = vr. 



Fig. m whose shaded region appears for sin^ 26*13 < 10"^ In Figs. [5]^c) and (d), the most 
stringent bound is obtained from fi 67. Note that sin^ 26*^^^^ in Figs. [5](c) and (d) become 
larger than 0.14 of the CHOOZ bound |l^ for mi > 0.008 eV, and then we can not have 
rrie^ = for the case. 

Figure [6] shows which process gives the most stringent lower bound on v/\mH±± in a 
space of Majorana phases by keeping = for the normal mass ordering. Green circles, 
blue crosses, and magenta triangles show regions where the strongest bound comes from 
/i — s> 67, r fiee, and the muonium conversion, respectively. It is impossible to achieve 
^efi = outside of the regions because of unacceptably large 9^^^, and then it becomes 
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(a) m, = O.OSeV, normal, m^ttlm^* = 1, rrig =0 
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(b) m-, = 0.2eV, normal, m^+->-/m^+ = 1 , mg^=0 
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FIG. 6: Green circles, blue crosses, and magenta triangles show regions where the strongest lower 
bound on v/\'mu±± for the case of = in the normal mass ordering comes from /i 67, 
r /xee, and the muonium conversion, respectively. The values of 9^^^ and d^^'^ are functions of 
Majorana phases (See Appendix). We can not have = outside of these regions. The regions 
are symmetric under a transformation (^^1,(^2) {—^i-; —^2)- (a) with mi = 0.05eV. (b) with 
rui = 0.2 eV. With rui = 0, the most stringent bound is given by /x ^ 67 for all values of (pi and 

^2- 



the situations discussed in the previous subsection for BR(/x eee) 7^ 0. Figure E](a) is for 
mi = O.OSeV and (b) is for mi = 0.2 eV. For mi = 0, it is possible to have = for 
any values of Majorana phases, and the most stringent bound is always given by // — ^ 67 
as v/^mn^^ ^ 200eV-GeV. Values of lower bounds on v/^mn^^ for Fig. [6(a) from fi ^ e^y 
and r /2ee vary in 150 350eV-GeV and 150 400eV-GeV, respectively. We have 
VAmH±± > 300eV-GeV, > 300 ^ 1500eV-GeV, and > 300eV-GeV in Fig. EKb) from 
fi — ^ 67, r — ^ /xee, and the muonium conversion, respectively. The regions are symmetric 
under a transformation (9:^1, 9^2) {—^1,-^2) because of |m^£/| = We see that 

991 ~ is preferred to keep s^^*^ small for mi 7^ and we can confirm s^^^ oc Am^i for 
9^1 = with eq. f|32)) in Appendix. It can be found also with eq. fl32l) for ipi,(p2 ^ 1 and 
Amf. = that 

ip2 ^ ^ '^^^'^ ifi = 0.68 ifi (red dotted line in Fig. EJ^b)) (29) 

is preferred to have a small s^^'^ . Bounds from r fiee and the muonium conversion can 
be the most stringent one only for ifi < O.Itt. Majorana phases are almost restricted as 
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FIG. 7: Lower bounds on vattt-h^^ foi' cases of nie^ = in the inverted mass ordering. The 
value of 6^^^ varies to keep niefj. = (See Appendix), (a) for (931,(^2) = (0,0), 5™^'^ = 0. (b) for 
{^1,^2) = {0,71), 6^^' = 0. 
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FIG. 8: The (^i-dependence of lower bounds on v/\mu±± for cases of rfig^ = with ms = in the 
inverted mass ordering. 6''^^'^ and (5™sc fQj- riif^^ = are functions of ipi (See Appendix). 



if 2 ^ 0.68 for the case of a strong constraint from the muonium conversion. It is shown 
that yU — > 67 can be the most stringent bound even for mi > 0.008 eV (cf. Fig. E]) because 
of large s^^'^ . At the border to the white region in Fig. El we have sin^ 2e'lf'' = 0.14. 

Similarly to Fig. [5l lower bounds on VATnH±± are shown in Fig. [7] for the inverted mass 
ordering. Note that ipi = it can not give = for the mass ordering because s^^'^ becomes 
too large. Very roughly speaking, the results in Figs. [Tl^a) and (b) are the same as those 
in Figs. M^a) and (b), respectively. A difference is that r fiee can be the most stringent 
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5 = 92 = 0, normal. m^*^/m^* = 1 , mgg=0 

2o C.L. <7^ with T data 

fi->ey 
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sin^2e,, 
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FIG. 9: Lower bounds on v^mfj±± for cases of iriee = with 6 = Lp2 = in the normal mass 
ordering. Values of m™^'^ and (p^^^ are given by the condition mge = depending on ^13 (See 
Appendix). 



bound at = while it is not the case for mi = in the normal mass ordering. Figure [8] 
shows the (/:>i-dependence for rrie^ = with = where the v?2-dependence vanishes. For 
ipi > O.Itt, sin^2^5^^^ becomes larger than 0.14. We see that /i ^ 67 can give the most 
stringent bound even for the inverted mass ordering. For 7723 = 0.2 eV, Fig. [6](b) is almost 
applicable to see which process gives the most stringent bound. 

Figure [9] is the result for the case of rriee = which is possible only in the normal mass 
ordering. Formulae of the magic values m™^'^ and yj™^'^ for = are shown in Appendix. 
The bound from /i — s> 67 is the most stringent one except for sin^2^i3 < 0.04 where the 
bound from r — > fififi becomes stronger than that. This is also the case for different values 
of 6 and ip2- 



III. CONCLUSION 



In the HTM, it is impossible to have a contribution to the muon anomalous MDM with 
a plus sign. Therefore, the HTM is qualitatively disfavored by positive Aa^Jr] (Aa^[e"''e~]) 
at 1.8cr (3.7(t). If we deal with 2a bounds to avoid the disagreement with Aa^Jr], v/^mH±± 
must be greater than 10^ eV-GeV in most of parameter space of the HTM in order to satisfy 
a strong constraint on BR(/i eee). We found that the bound on vi^mfj±± from r Jiee 
becomes more stringent than that from /i — > eee in a region of nii > 0.06 eV, ipi < 0.0027r, 
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O.Svr < ipi < I.Stt, and sin^ 2^13 < 10~^. The bound from /i eee can be evaded in cases 
of rrie^ = and niee = 6|. In the case of = in the normal mass ordering, the 
strongest bound is given by /i — ^ 67 for mi < 0.01 eV or ifi > O.lvr and by r decays (mainly 
r /2ee) for mi > 0.01 eV with Lpi < O.lvr. On the other hand, r fiee gives the strongest 
bound in the normal mass ordering except for m^, ~ with ipi ~ O.Itt where fi ^ e'j gives the 
bound. For both of the mass orderings with mg^ = 0, the muonium conversion gives the most 
stringent bound if Majorana phases satisfy (p2 — O.GSipi for (pi,(p2 -C 1 and mi > 0.1 eV. 
In the case of mge = 0, the strongest bound is obtained from /i — > 67 except for the case of 
sin^ 2^13 < 0.04 where the bound is given by r — /i/i/i. By looking over all cases, we see that 
V/\jnfj±± > 150eV-GeV should be satisfied in the HTM. If mH±± is measured, the bound 
can be the lower bound on v/^ though there remains a possibility of v/^ = for which we 
can not use the correlation of hui with \^V/\hii' = mw. 
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APPENDIX 

We obtained formulae of the magic values of ^13 and 5, which give mg^ = 0, as 
cos^^g^ - C12S12C23C 

~ mgc (2 42-2 / 2 i2 ^2^' ^ ^ 

■523'5i3 (mg — Si2m2 Sm ipi — (512^2 COS 991 + C^2"^lJ / 

sin<5-g= = C12S12C23D .g^x 

-323-513 |m3 — 312^2 V^l ~ (-5i2"^2 COS^^i + Ci2miJ^| 



S 



mgc 



C12S 12C23 Vc^TlT' 



S23 |m| — s\2m\ sin^ (fi — (sf2"^2 cos (pi + c\2^iY I 



C = —Ci2mi + Si2m2 + mim2 cos 2^12 cos v^i 



(32) 



- mim3COSV92 + m2m3 cos(v9i - V92), (33) 
D = — mim2 sin <y5i + mim3 sin (y92 + m2m3 sin((y9i — ^92), (34) 
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where we define s^^'^ = sin 9^^'^. Note that s^^^*^ 7^ because it requires m2 = mi. Note also 
that s^^'^ is not always acceptable; For example, = vr in the inverted mass ordering gives 



-.mgc 
'13 



> 1. These results are consistent with s^^'^ and 6'^^'^ used in 



Although mixing matrix in 



14| for (y92 = or TV. 



6[ is defined as = U^ii^i in stead of Vi = UuVi used in 
this article, there is no change in formulae of magic values because the difference appears 
just as the simultaneous flip of signs of all phases. 

On the other hand, = can be achieved only in the normal mass ordering. The 



magic values of (pi and mi for mge = are given as functions of S13 and — 25 6| by 



mgc 

sm 



,/{mf^y + Am: 



31 



^i2\/ (''^T^'^y + Am 



tl^sm{ip2 — 25), cosip^^'^ < 0, 



(35) 



m 



mgc\ 2 
1 ) 



C0S2 2^12 - 2 {Sf2 + Cf2 COS 2(v?2 - 25)) tfg + tfg 

S12 cos26'i2Am2i + {siaAmgi + (si2 + C12 cos 2(^92 — 25)) Aml^^ t 



X 



4 

13 



A ^ (4A 
B 



m. 



21 



- Amg^t^g - 2012^13 cos((/?2 -25)^ A + Bt\^ , (36) 

cos 26'i2Am3;^) s^2^"^2i5 (37) 

(38) 



{(■5^2 - C12 sin^((y92 - 25))Am3i - s^2Am2i| Amg^, 
where we define tn = S13/C13. 
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